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Abstract 

We study general constraints on spontaneous R-symmetry breaking models coming from the 
cosmological effects of the pseudo Nambu-Goldstone bosons, R-axions. They are substantially pro- 
duced in the early Universe and may cause several cosmological problems. We focus on relatively 
long-lived R-axions and find that in a wide range of parameter space, models are severely con- 
strained. In particular, R-axions with mass less than 1 MeV are generally ruled out for relatively 
high reheating temperature, Tr > 10 GeV. 
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I. INTRODUCTION 



Supersymmetry (SUSY) has been considered to be the strongest candidate of the physics 
beyond the standard model (BSM). Although the recent data from the Large Hadron Collider 
(LHC) has not shown any evidence for SUSY but discovered a Standard Model Higgs-like 
particle with a mass of around 125 GeV it still remains a strong candidate of BSM 
because it suggests the gauge coupling unification, it guarantees proton stability sufficiently, 
and it provides a reasonable dark matter candidate. Moreover, in string theories, which are 
the most powerful candidates of the quantum theory of gravity, it plays a crucial role for 
consistency and must be broken at a scale between the electroweak scale and the Planck 
scale. _Theref ore, it is important to investigate SUSY-breaking models in the light of LHC 
data 

R-symmetry, which is a specific symmetry of supersymmetric models, is a key ingredient 
for SUSY breaking and its application to model building. Recent drastic progresses on SUSY 
breaking by exploiting a metastable state (see js-^l for reviews and references therein ) give s 
u= a bette. unde.tanding otthe .ole of the R-syo,™t.y in realistic .nodel building M 
Nelson-Seiberg's argument [20[ beautifully demonstrates a connection between metastability 
and R-symmetry in the context of generalized Wess-Zumino models with a generic super- 
potential. If R-symmetry is preserved, there is no SUSY vacuum in a finite distance in 
field space. On the other hand, if a gaugino mass has Majorana mass, R-symmetry has 
to be broken to generate the gaugino mass. Thus, there is a tension between stability of 
vacuum and generating gaugino mass. A simple solution to this problem is to introduce an 
approximate R-symmetry. 

One of the interesting ways to break R-symmetry is spontaneous breaking. In Ref. 2l|, 



D. Shih revealed a quite fascinating condition for spontaneous R-symmetry breaking in the 
context of generalized O'Raifeartaigh models: For R-symmetry breaking, there must be a 
ield with R-charge different from or 2. Such models were applied to gauge mediation 
22[ and some classes of the models successfully generated large gaugino masses. According 
to the general argument by Komargodski and Shih 23|, large gaugino mass is related to a 
tachyonic direction at a point in pseudo moduli space toward the messenger direction. In 
the R-symmetric model, such tachyonic direction exists at the origin of the pseudo moduli 
space. 
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When the spontaneous breaking of ?7(1)_r symmetry occurs, cosmic R-strings are formed 



by the Kibble-Zurek mechanism 



24 



251 ] ■ Plugging the structure of the pseudo-moduh space 



mentioned above and R-string forming, we will meet a quite dangerous possibility. It is 
known as a "roll-over" process of vacuum through inhomogeneous energy distribution by an 



impurity such as a cosmic string 26|, l27(|. In the core of the R-string, the system can easily 



slide down to the lower vacuum via the tachyonic direction at the origin and form a sort of 
"R-tube" in which the core sits in the lower energy vacuum. Thus, if the tube is unstable, by 



rapid expansion o 
discussed in Ref. 



Ref. 
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the radius, the universe can be filled by the unwanted SUSY vacuum. As 



28| this gives a constrain for model building. However, as emphasized in 



30|, when a D-teim contribution is not negligible, it can lift the tachyonic direction 



and stabilize the pseudo-moduli space. In such models, the roll-over process does not occur. 
Also, when the amplitude of (tachyonic) messenger mass at the origin is sufficiently smaller 
than that of R-symmetry breaking field, the vacuum selection is successfully realized. As 
we will see, R-strings are unstable due to the explicit R-symmetry breaking term in the 
superpotential and hence the roll-over process can be circumvented if the life-time due to 
the explicit R-symmetry breaking is shorter than that for the roll-over process. In this paper, 
we assume such an early stage scenario and study general cosmological constraints for the 
models. In this sense, the results shown in the present paper is complementary to the ones 
studied in Ref. 



In spontaneous R-symmetry breaking models, there exists a pseudo Nambu-Goldstone 
boson, called R-axion, as well as the modulus field called R-saxion. They are copiously pro- 
duced in the early Universe from scattering of thermal plasma, coherent oscillation, R-string 
decay and so on, and may cause other cosmological problems. Therefore, model parame- 
ters on spontaneous R-symmetry breaking model can be constrained from such cosmological 
considerations. Note that, unlike the QCD-axion, R-axions receive relatively heavy mass 
from gravitational coupling with explicit R-symmetry breaking constant term in the super- 
potential and its lifetime can be much shorter than the cosmic age. Thus, we can impose not 
only constraints from the R-axion overclosure problem but also that from R-axion decay. 
In this paper, we investigate their cosmological constraints focusing on relatively long-lived 
parameter range. We show that the model parameter space is severely constrained and many 
parameter space of R-axion is ruled out from the cosmological consideration. 

This paper is organized as follows. In section [TTl we explain the general feature of spon- 
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taneous R-symmetry breaking models. In section IIIIl we evaluate the R-axion abundance 
produced in the early Universe. Here we assume that cosmic R-string is produced in some 
earlier epoch. We list the cosmological effects induced by R-axions in section IIV[ We also 
evaluate the constraint on the parameter space from these effects. Section |V] is devoted to 
conclusion and discussion. 



II. SPONTANEOUS R-SYMMETRY BREAKING MODEL 

In spontaneous R-symmetry breaking models, the SUSY-breaking field with a finite R- 
charge acquires nonvanishing vacuum expectation value. The phase of the SUSY-breaking 
field is almost massless and identified as the Nambu-Goldstone boson. It acquires a small 
mass from explicit R-symmetry breaking term in the superpotential and called R-axions. 
In order to see its cosmological consequeces, we should first investigate their properties and 
interactions. Here, we review a simple but general R-symmetry breaking model focusing on 
R-axions and read off their interactions with several modes. 



A. R-symmetry breaking model 

Let us consider a simple effective superpotential for the R-charged SUSY-breaking field, 
X, integrating out the messenger fields, 

W,s = AlsX + Wo. (1) 

Here Aofj gives the nonvanishing F-term for the SUSY-breaking field and R-symmetry break- 
ing constant Wq is introduced for the cosmological constant to vanish. Note that from the 
flat Universe condition, they are related as A^g = SW^/M^^ with Mpi being the reduced 
Planck mass. Assuming a noncanonical Kahler potential, X can be destabilized at the 



origin 



Here we consider the effective potential for X, 
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V{X) = ^{\X\^-f!f + '^f.X + h.c. 

= ^ix' - VI? + ^/aXcos(a/v^/„), (2) 

where we have defined X = (x/A/2)e*'^/^^''. The second term that breaks U{1)r symmetry 
comes from the R-symmetry breaking constant term in the superpotential that couples to 



X field through the Planck suppressed interaction in supergravity. The R-axion mass is 
related to the parameters in the potential as 

" ~lMl ~ Ja ' 

where 7713/2 = Wq/M^^ is the gravitino mass. 

Let us investigate the model further. We here expand X around X = fa as follows, 

X = ^±^exp(za/v^/.), (4) 

so that the fields a and s have canonical kinetic terms. Here the phase part a and the radius 
part s are identified as R-axion and R-saxion, respectively. Note that the mass of R-saxion 
is related to the R-symmetry breaking scale as 



"Z^fa^ \l^rna. (5) 

Ja 



In the last equality, we assumed that the Kahler metric is given by 



9xx-^^-p\^\' + T^\X\^ (6) 

J a J a 



with A being a numerical constant of order of the unity, and A is related to the model 
parameters as A ^ m^Mpi//^. The fermionic partner of X field, "R-axino," is the goldstino 
for the SUSY-breaking and absorbed in the gravitino. Thus, we just have to consider 
cosmology of gravitinos instead of R-axinos. 



B. Interactions of R-axions 

We now investigate interactions of the R-axion with several modes as well as its cross 
sections and decay rates. As we will see, they are useful for the cosmological constraints on 
R-axion abundance. 

First of all, the R-axion to R-saxion interaction can be read off from the kinetic term of 

X, 

|5.^p4(l + ^)'(9,a)l (7) 
From this interaction, we can evaluate the decay rate of R-saxion to 2 R-axions as 



We can assign R-charges to the supersymmetric Standard Model fields such that the 
R-symmetry is consistent with all of the interactions. After SUSY and R-symmetry are 
broken, the R-axion appears in the gaugino mass terms as well as the so-called B-term and 
A-terms. In addition, the R-axion couplings with the gauge bosons appear through the 
anomaly coupling terms. That is, the coupling between the R-axion and the photon is given 
by 



lem 



aF^,F^\ (9) 



where F^.^ is the field strength tensor of U{\)era and Cem is the anomaly coefficient, i.e. 
Tr ?7(l)i{f/(l)g^, which is mo del- dependent. Then, the decay width of the R-axion into two 
photons is given by 

167r V 47r y V /a 



■p/ V o \ ^ em ( 9em\^ f IT^a i 

T{a 27) ~ -— — [-T ] 



Similarly, the R-axion coupling with the gluon is given by 

T^aG,.G'^^ (11) 

where G^^ is the SU(3) field strength tensor and Cg is the anomaly coefficient, i.e. 
Tr U{l)jiSU{3)^. R-axion decay into gluons are associated with QCD jet production and 
hence R-axion must be heavier than at least the proton mass for this decay channel, which 
is beyond our interest. Thus, we do not consider it here. 

The interactions of the R-axion with the Higgs fields appear through the B-term. Then, 
the R-axion and the Higgs fields mix each other in their mass terms (see for its detail 
Appendix lAl ) . The eigenstate corresponding to the low-energy R-axion a includes the axial 



parts of the up and down-sector Higgs fields, and [31 1 



a ~ Kfl + Kr cos^ (3 sin /3 + Kr sin^ (3 cos (3 (12) 

where r = v/ (y^/a), v = 246 GeV, k = (1 + r^ sin^ 2/3)"-'^/^. Note that a denotes the R-axion 
at high energy beyond the electroweak symmetry breaking. Since the coefficients of ^u.d are 
very small, it is found that a ~ a. Hereafter, we denote the low-energy R-axion by a instead 
of a. However, because of this mixing, the R-axion can couple with the quarks and leptons 
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through their Yukawa couphngs. That is, the couphngs of the R-axion with the up-type 
quarks, the down- type quarks and the charged leptons, Au, A^ and A^, are given by 

Am = ii/ui^r cos^ f3 sin (3 = i—^K cos^ (3, 

J a 

= ^j/d^'^sin^ /3cos/3 = z— ^Ksin^ /3, (13) 

Ja 

\ • 'In n • • 2 n 

A^ = lygKr sm pcosp = t—Ksm p, 

J a 

respectively, where yf and rrif with / = u,d,£ are their Yukawa couplings and masses. 
Through these couphngs, the R-axion can decay to a pair of the SM fermions, if rria > 2771 f. 
Its decay width is given by 

r(a ^ //) = (1 - 4mJ/m^) . (14) 

For example, the decay rate into the electron pair is given by 



1.0 X 10-3iGeV X sin^/3 (^) (^^^)' {l - 4ml/mlY^\ (15) 

The decay rate into the fi pair is enhanced by its mass as r(a — /i'''/i~) = {fn^/meY x r(a — )■ 
e^e~), but such a decay occurs for nia > 2m^. 

Similarly, we can compute the couphngs between the R-axion and the neutrinos. For the 
neutrinos, we consider the Weinberg operator in the superpotential, yy{LHuY / Mr, instead 
of the Yukawa couplings terms. Then, similar to the above, the couphng of the R-axion with 
neutrinos is given by 

= i—^Kcos^ (3. (16) 

Ja 

Thus, the decay rate of the R-axion into the neutrino pair is suppressed because it is pro- 
portional to the neutrino mass squared, i.e. 

/ \ 2 



T(a-^uiy) = { — ] cot^ (3 xT(a^ee). (17) 

Therefore, the branching ratio of R-axions into pair of neutrinos are small enough even in 
the case where the decay channel into electron is closed, rria < MeV. 



III. R-AXION PRODUCTION IN THE EARLY UNIVERSE 



Let us consider the cosmology of the spontaneous R-symmetry breaking model focusing 
on the R-axion production and evaluate the R-axion abundance. We consider the case 
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where U{1)r is restored due to some additional mass terms such as the Hubble induced 
mass or thermal mass in the early Universe^. After some epoch, X field is destabilized as 
the additional mass term decreases and acquires vacuum expectation value x = fa- Since 
the approximate f/(l)/? symmetry breaks spontaneously at that time, (unstable) cosmic 
strings are formed by the Kibble- Zurek mechanism. The long cosmic strings in a Hubble 
volume intersect each other and generates closed string loops^. These closed string loops 
shrink with emitting R-axions. As a consequence, the cosmic string network enters the 
scaling regime. As the Hubble parameter decreases, the explicitly R-symmetry breaking 
term in the potential becomes no longer irrelevant to the dynamics of the system and the 
string network turns to the string-wall system where domain walls are attached to cosmic 



strings 
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33|. The string- wall networks are unstable and annihilate when the domain wall 



tension becomes comparable to that of cosmic strings. The energy stored in the string-wall 
system turns to R-axion particles. The fate of R-axions produced from the cosmic string 
loops and the string-wall system as well as the scattering of thermal plasma and the vacuum 
misalignment is determined by the lifetime of R-axions, which, then, constrain the model 
parameters of spontaneous R-symmetry breaking models^. In the following, we estimate the 
R-axion abundance from each source. We will examine the cosmological constraints in Sec. 



A. R-axion production from vacuum misalignment 



First we evaluate the energy density of the coherent oscillation of the R-axion field 34 1. 
After the spontaneous R-symmetry breaking phase transition, the R-axion field acquires 
some initial value, a^, and keeps its position after a while due to large Hubble friction. 



^ We assume that the SUSY-breaking vacuum is selected by some mechanism. Note that if the amphtude of 
the mass of X at the origin is larger than that of messenger fields, the SUSY-breaking vacuum is naturally 
selected. However, this issue is beyond the scope of this study and hence we do not impose any constraints 
on the model parameters from the vacuum selection. 
^ Here we do not consider the effect of the existence of SUSY vacuum on the cosmic string structure. This 
issue will be studied elsewhere 



28|. 



R-axions are also produced from R-saxion decay. However, as shown in Appendix |B1 the abundance of 
such R-axions are generally subdominant and hence we do not consider it here. 
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When the Hubble parameter decreases to the R-axion mass, 

Hitosc) = rua, (18) 

the R-axion field starts to oscillate. Here the subscription "osc" indicates that the parameter 
or variable is evaluated at the onset of the R-axion oscillation. The energy density of the 
oscillating R-axion pa^o is given by 

paAtosc) = \mla\. (19) 

If the R-symmetry is broken after infiation, the initial value of R-axion distributes randomly 
from —\plT\ja to \pl'Kfa siucc the correlation length of R-axion becomes much shorter than 
the Hubble length at the onset of the R-axion oscillation. Therefore, we estimate the mean 
value of as 

(a,2) = — =— / a\da. = — ^. (20) 

Since the energy density of R-axion oscillation decreases as a~^, the quantity pa/s is con- 
served as long as there are no entropy production, where s is the entropy density. Therefore, 
we characterize the axion abundance by this quantity as 

15 mlfa 



for Hcsc < H 



R 



f f' ^ Th. for > H„ 



where g^, and g^:s are (effective) relativistic degrees of freedom for energy density and entropy, 
respectively, and the subscript "i?" represents that the parameter or variable is evaluated 
at reheating. Note that Tqsc is given by 

^ (to) ' (iSv) " ■ 

Here we assume that the scale factor increases like matter dominated era during infiaton 
oscillation dominated era and take into account the dilution until the infiaton decay or 
reheating when i^osc > Hji. 



B. R-axion production from global cosmic strings 



Next we evaluate the energy density of R-axions radiated from the cosmic string loops 
35l . |36| following the discussion in Appendix B of Ref. 37|]. When the R-string network 
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enters the scaling regime, the energy density of the long R-strings are estimated as 

M«) = l^/^ln (^) . (23) 

\ '^string / 



Here the scaling parameter ^ ~ 0.9 37|, |38|] represents the mean number of strings in a 



Hubble volume and (istring — ^ ^^'^fa ^ represents the core width of R-string. Note that the 
line energy density or the tension of R-string is given by 



^ 2vr/^n f y^^y (24) 

\ 'J^string / 

Assuming all the energy loss of long R-strings is converted into R-axion particles through 
the string loops, we obtain the evolution equations 

dpoo{t) 



dt 

dpa,str{t) 



dt 

where the energy emission rate from the string loops 



-2Hp^{t)-r,^{t), (25) 

-AHpa,stM+^cmit), (26) 







'^string / 
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V '^string 



^ ^ .2 In -j-^ - 1 ' for RD 

re.(t) = ^x<| )^ ^^/,7.4x \ (27) 



Here we assume that R-axion particles released from cosmic string loops are relativistic. 
Since the mean comoving momentum of radiated R-axion can be evaluated as 

ka,str{t) _ 2ne 



(28) 



R{t) t 

with the constant e ^ 0.25 [36|, |37[, we can estimate the number density of radiated R-axions 



as 



R{ty Jt, KMt') 

'^^fa V J V V ^string/ / /onN 



In ( 1 - 3 1 , for t > 



3 \ \ C^string / / 

Here is the time when the R-string network enters the scaling regime. 

When the Hubble parameter becomes comparable to the R-axion mass and R-symmetry 
breaking mass term becomes no longer irrelevant, t = tosc string- wall system forms and 
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R-axion emission from R-string loops stops. We can evaluate the resultant number density 
of R-axions from the R-string loops as 

frr^ I^P^f^^ ^ )-3),for H,,,<Hr 



The radiated R-axions become nonrelativistic after some epoch. Therefore, we can approxi- 
mate the R-axion energy density as Pa,str = ''^a'^a.str and the R-axion energy-to-entropy ratio 
as 

90 e </aV, f 1 



Pa.str _ J TT^ g*s{Tosc)e T^^^ V^'^ \2maV^4tring ^ ' Hose < Hr 



4fi'*.(TR)e V^pi/ V \3mav^4tring/ 2, 
C. R-axion production from string- wall system 

Let us evaluate the energy density of R-axions from the string-wall system annihilation 
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33| . At t ~ tosc; the explicitly R-symmetry breaking term in the potential becomes 
no longer irrelevant, and string- wall system forms. The surface mass density of domain walls 
are estimated as [39| 

o-„aii = IQrriafa- (32) 
When the tension of domain walls dominates that of strings. 



c^waii = =1 ^ t In ( ^ ) = '-m-^\ (33) 



we assume that 



the string-wall system annihilates. As following the discussion in Ref. 
the energy stored in the string-wall system released to R-axion particles. Thus, we evaluate 
the number density of R-axions as 

Pwall osc ) + Pooitosc) ( Rit 



'^a,sw (^) 



0"wall ^ Pstring(^osc) | / -R(^osc) | (34) 



Ot-^Tfld Y ^osc ^osc / V Ri.^^ 



where Ua = a^rna is the average energy of radiated axions and A = p„aii^/^waii — 0.5 (33| is 



the area parameter of domain walls. The radiated R-axions become eventually nonrelativistic 
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and hence we can evaluate the energy-to-entropy ratio as 

180 m2/2 



\ 2 (T \ T3 \ 7^ Tfl ) ) ' for ^osc < Hr 

Pa,sw fT^a^afivf J ^ fi'*s osc J C^io osc \ \ ^'^aV S '^string , 



^^"^V.f24^+^einf , i )).for H^..> Hr 



^string 



(35) 



string j — 



Noting that the logarithmic factor is evaluated as \n{l /ma^/^d, 
\n{^/X/^{fa/ma)) = 30 for ~ 10^° GeV and m„ ~ 1 MeV, hereafter we approxi- 
mate the R-axion abundance from R-axion dynamics, i.e., the coherent oscillation, the 
decay of cosmic string loops, and the decay of the string-wall system, 

n I ^l^' H^sc < Hr 

f a,dyn ra,o i pa,str i ya,sw ) 1 Qg(, 



9,4xlO-»GeV/f,(^)"V^) . for 



^ ilo4v) (lo^) • > ^« 

where Ki ^ 0(1) and K2 — 0(10) are numerical parameters. 



(36) 



D. R-axion production from thermal bath 

We have estimated the abundance of R-axions generated from their dynamics. We should 
also take into account that generated from other sources. Here we evaluate the R-axion 
abundance from thermal bath. The R-axion abundance from R-saxion decay is discussed in 
Appendix [B] and is generally negligible. 

R-axions are produced in the thermal plasma from (mainly) gluon scattering, gg — )■ ag. 
Since the gluon-axion interaction comes from the anomaly term, 

with Cg being the model dependent anomalous coefficient and gs being the strong gauge 



coupling, the R-axion abundance is calculated as [40 



Pa,ih 



^■°-°-'«--n^(l^)(™)'(l0^)- (38) 
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Note that R-axions are thermalized once if the reheating temperature is high enough, 

Tn>Tn^ 10^GeV,;^C- (jj^) ' , (39) 

where is R-axion decouphng temperature. In this case, the R-axion abundance is evalu- 
ated as 

^ ^ 2.6 X lO-^GeV ( . (40) 
s VlMeV/ ^ ^ 

Note that R-axion is produced thermally only if Tr > is satisfied. 

As a result, the total R-axion abundance in the early Universe is evaluated by the sum 
of these contributions and given by 

Pa _ Pa, dyn ^ Pa,th (41) 
S S S ' 

In Fig. [H we show the theoretical predictions for the R-axion to entropy ratio with 
fa = lO^GeV, lO^GeV, 10^°GeV, and lO^^Qgy. Here the solid lines represent contribution 
from the thermal production (Eqs. ( |38l) and (l40l) ) and dashed ones represent the R-axion 
dynamics (Eq. fl5^ ) with Ki = 1, K2 = 20, respectively. Black, blue, green and red lines 
correspond to Tr = lO^^GeV, IGeV, lO'^GeV, and lO^GeV, respectively. In the case of 
Tr > Td, R-axion abundance from thermal production is independent of Tr. We can see 
that the contribution from the thermal production generally dominates for < lO^^GeV. 
Nonthermal contribution dominates for smaller and larger /„. 

IV. COSMOLOGICAL CONSTRAINTS FROM R-AXION 

Now we consider the generic constraints of the R-symmetry breaking model from cosmol- 
ogy. One may think that the model with long-lived R-axions is safe if they never dominate 
the energy density of the Universe or R-axions are responsible for the dark matter in the 
present Universe. However, even if they are subdominant component of the Universe, their 
(partial) decay is constrained by several cosmic/astrophysical observations depending on 



their abundance j43|. Since we have evaluated the R-axion abundance and its lifetime, we 
can constrain the model from various observations. As we will see, strong constraints for 
the model parameters are imposed. 
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FIG. 1: Theoretical predictions for the R-axion to entropy ratio with fa = lO^GeV, lO^GeV, 
and lO^^GeV. The sohd hnes represent contribution from the thermal production 
(Eqs. ([38]) and (|30]) ) and the dashed ones represent the R-axion dynamics (Eq. (j36|) ) with Ki = 1, 
K2 = 20. Black, blue, green and red lines correspond to Tji = lO^^GeV, IGeV, lO^GeV, and 
lO^GeV, respectively. 

A. Constraints constraints on axion abundance 

Let us see the various constraints of R-axion abundance from cosmology and astrophysical 
observations. 



1. Big Bang Nucleosynthesis 



The R-axion decay into photon or electron (radiative decay) after the Big Bang Nu- 
cleosynthesis fBBN) epoch may break the light elements and the R-axion abundance is 



constrained 



44| . The radiative decay of R-axion causes photo-dissociation process of light 



elements and changes the light elements abundance. We can read off the constraint on the 
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R-axion abundance at its decay from Ref. 4J] as 



n ( lO^^GeV ( , for lO^s < < lO^s 
B,^<{ V10%y (42) 

^ ^ ( 10-i^GeV, for lO^s < < lO^^s 

where B,- is the radiative branching ratio^. Note that this effect is neghgible if the energy of 
the injected photons is so small that they cannot destroy the light elements. Thus, we here 
impose a condition for this constraint to be effective, 

> 4.5MeV, (43) 

which corresponds to the threshold energy for the deuteron destruction process, -D + 7 — )■ 
n + p. 



2. Cosmic microwave background distortion 



The radiative decay of R-axion before the recombination may distort the blackbody 
spectrum of CMB. After the double- Compton scattering freezes out at t ~ 10^ s, energy 
injections generate nonzero chemical potential fi of the CMB spectrum, which imposes the 
constraint from the blackbody spectrum distortion of CMB. Energy injections after t ^ 10^ 
s, when the Compton scattering is no longer in thermal equilibrium, thermalize electron, 
which causes the Sunyaev-Zel'dovich (SZ) effect. Since the SZ effect is constrained by the 
Compton y-parameter, we can impose a constraint on the R-axion abundance. 



The COBE FIRAS measurement 



451] constrains the CMB distortion as 



|/i|<9xl0-^ y<1.2xl0-^ 



(44) 



Since the injected energy is related to these parameters as [46|, |47| 



^ ~ 0.714/i, for lO'^s <Ta< lO^s 
^ ~ 4y, for lO^s < < lO^^s 

P7 



(45) 
(46) 



^ If the R-axioii mass is heavy enough, ma > 2 GeV, the hadronic decay channel opens. In this case, more 
stringent constraints are imposed |44l |. 
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the constraints on the R-axion abundance is given by 

< lO-i^GeV (^^y^' for lO'^s < r„ < lO^^g (47) 

depending on its hfe time. Note that /i and y-parameters impose almost the same constraint 
on the R-axion abundance at its decay. 



3. Diffuse X-ray and j-ray background 



The R-axion decay to photons after recombination, t > lO^"' s, may be constrained 
from the diffuse X-ray and 7-ray background observation. Photons with energy IkeV < 

< ITeV rarely scatter with the CMB photons and intergalactic medium. Therefore, the 
photons produced from the R-axion decay in the "transparency window" 48[, 



r > 



10^9 



rria y 



keV/ 



4 X 10^^ s. 



nir, 



100 MeV 



for 1 keV < < 100 keV 

for 100 keV < < 2.5 MeV 

for 2.5 MeV < < 100 MeV 

for 100 MeV < m < 10 GeV 



propagate through the Universe and can be detected as diffuse background. 
The flux of the extragalactic diffuse photons is roughly given by 

-0.4 



2 X 



E 



(^;/30keV)0-3 + {E/30keVy- 
MeVy 



9' 



5.0 X 10- 
1.7 X 10"^ 



-1.1 



1.45 X 10- 



VlOOMeV J 

\ —1.4 



( 



E 



V lOOMeV 



0.25keV <E < lOkeV 
lOkeV <E < SOOkeV 

SOOkeV <E < 30MeV 

30MeV <E < lOOMeV 

lOOMeV <E < lOOGeV 



Here we applied the observational results of ASCA j49| for 0.25-10 keV, HEAO [501 fo'^ 25 



(48) 



(49) 



keV- 800 keV, COMPTEL |5l| for 800 keV-30 MeV, EGRET [52] for 30 - 100 MeV, and 
Fermi 53| for 100 MeV-100 GeV. Note that we have taken into account the resolved source 



of diffuse X-ray background 5^, l55|] and used the fitting formula derived in Ref. 56 1. 
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The flux of photons produced from the R-axion decay can be approximated as 

na.o 



F,{E) 



~ i?^ X 



2nTaHo ' 

3?T.a,dec So 



for Ta > to 
for Ta < to 



(50) 



Air Sde 

where the subscriptions "0" and "dec" indicate that the parameter or variable is evaluated 
at the present and the R-axion decay time, respectively, and is the branching ratio to 
photons. Note that the energy of photons should be evaluated at E = ma/2 for Ta > to and 
E = (3i/oTaV^m/2)^/'^(ma/2) for Ta < to, taking into account of the redshift of the photons. 
Then, the abundance of the R-axions are constrained from the constraint F^{E) < F°^^ as^ 



B — < 

s 



2Ah X 10-i^GeV 



4.8 X IQ-^^GeV 




F°b«(m,/2) 



10 ^cm ^s ^ 




(51) 



,lMeV/ I 10-W-2s-2 
where h = iJo/(100 km sec~^Mpc~^) and Ho is the present Hubble parameter. 



4- Reionization 

The radiative decay of R-axion after recombination is also constrained from reionization. 
If the energy of injected photons is relatively small, they are redshifted and interact with 
intergalactic medium. Then, the intergalactic medium is partially ionized and the R-axion 
decay is regarded as an additional source of reionization. To be consistent with the obser- 
vation of the optical depth to the last scattering surface, the R-axion abundance should be 
small enough. Assuming that the one-third of the energy of photons produced from R-axion 
decay that leaves the transparency window is converted to the ionization of the in 
medium, the R-axion abundance can be constrained from the inequality in Ref. 

f 6.77 + 3.96275X + 0.25858x2 + 0.00445^3, -17 < x < -13 
logioC< (52) 
y -24.75 -X, x<-17 

^ Most of diffuse extragalactic X-ray and 7-ray background can be explained by astrophysical sources such as 
blazers. However, here we use the conservative constraint, though the flux of the unresolved extragalactic 
defuse background photon would be much more smaller when we assume some cosmological models of the 
evolution of galaxies. 
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where ^ 

C = firPa/pbaryonldec = 0.43 X lO^^GeV"^ (^^^ (53) 

and X = log^o(r/s~^) = — log^Q(ra/s). Here fib denotes the present density parameter of the 
baryonic matter. This constraint is complementary to that from the diffuse X-ray and 7-ray 
background. 



5. Dark matter abundance 



If the hfetime of R-axions is longer than the present time to, most of R-axions remain the 
present Universe and contribute to the dark matter of the Universe. Thus, we can constrain 
the R-axion abundance in order not to exceed that of the dark matter. In terms of the 
energy-to-entropy ratio, the R-axion abundance is constrained as jsSi 

^ < 4.7 X 10-i°GeV (^^) • (54) 
B. Constraints on model parameters 

Now we are ready to show cosmological constraints for spontaneous R-symmetry breaking 
models. In FigJJj we show the constraints on the model parameters, rria and fa coming from 
various conditions argued in the previous subsection. Each colored region is excluded and 
white region is allowed. As a reference, we also show lines of gravitino mass. Upper dotted 
lines and lower ones represent ^3/2 = IkeV, 7713/2 = leV, respectively. Here we focus on the 
region > 10^ GeV since smaller fa is forbidden from laboratory experiments such as rare 



decays of or 59 1 



For the higher reheating temperature, Tji > 10^ GeV, all the parameter space where R- 
axions decay at t > 10^(10^) sec is ruled out for > 4.5 MeV (m^ < 4.5 MeV) regardless 
of reheating temperature. This is because R-axions are inevitably produced so much that 
cannot pass any constraints discussed above. Short-lived R-axion is allowed because any 
entropy production before is not forbidden and there are no cosmological constraints. 

On the other hand, for the smaller reheating temperature, Tr < 10^ GeV, several param- 
eter space where R-axions decay later is allowed. This can be understood from Eqs. ( I36ll 
and (138|) . For larger /„, nonthermal production is dominant and the R-axion abundance is 
expressed as pa/ s oc /f, whereas thermal production, which depends on fa as pa oc f'"^, 
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TR=10"^GeV 



Jp-lGeV 




VMeV/ 





-3-2-10 1 2 

Logiol — ^ I 
VMeV/ 

TR=10^GeV 



Logiol — ^ I 
VMeV/ 




FIG. 2: Cosmological constraints on the naodel parameters, nia and fa with Tr = 10 "^GeV, IGeV, 
lO'^GeV, and lO^GeV. Each colored region is excluded and white region is allowed. Upper dotted 



lines and lower ones represent — l^eV and m 



Hl2 



\eV. 



dominates for smaller fa- Thus, the R-axion abundance takes its lower value at fa ~ 10^^ 
GeV. As a result, allowed parameter region appears at rria ~ 10 keV and 1-100 MeV for 
fa ~ 10^2 Tn ~ 10^2 - 1 GeV. 

One may regard that there is a parameter region where R-axions can be dark matter for 
smaller reheating temperature, ~ 10^^ GeV. However, in this parameter region, R-saxion 
mass is considerably small, < 1 MeV for a naive model discussed in Sec. [Tll For the 
point of view of vacuum selection and R-string stability 28|], this requires unacceptably small 
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messenger mass. Therefore, we conclude that an ingenious model building is necessary for 
R-axions to be dark matter of the present Universe. 

Thus far, we did not take into account constraints from R-axinos or gravitinos. Gravitinos 
are produced from gluino scattering in thermal plasma and their abundance is evaluated as 
6o|, 

^.,.x:0-GeVx(^)=(^)-'(^), (55) 

where rrig is the gaugino mass. Since gravitinos are stable in this case, gravitino abundance is 
constrained from dark matter abundance (Eq. fl54l) ). As a result, depending on the gaugino 
mass, another stringent constraint is imposed for model parameters in higher reheating 
temperature case Tr > 10 GeV: Smaller fa and would be forbidden. 

In summary, we have shown that spontaneous R-symmetry breaking models are severely 
constrained from cosmological considerations and generally long-lived R-axions are forbid- 
den. In order to avoid that, careful model building and smaller reheating temperature are 
required. 



V. DISCUSSION 



We have studied general cosmological constraints on spontaneous R-symmetry breaking 
models. We estimated the abundance of R-axion produced firstly via their dynamics such 
as coherent oscillation and decay of cosmic string/wall system, and secondly via thermal 
scattering process from gluon-axion interaction. It is interesting that R-axion production 
from R-string and wall systems are large enough and can be dominant in some parameter 
region. Basically the models were motivated by gauge mediation, gravitino as well as R- 
axion are relatively light. Therefore, R-axion tends to be long-lived. The conditions for the 
R-axion density coming from BBN, X-ray/7-ray background, reionization and overclosure 
severely constrain the scale of R-symmetry breaking. As a result, smaller R-symmetry 
breaking scale and SUSY-breaking scale are disfavored from cosmological constraints. In 
the point of view of gauge mediation, this result weakens its motivation, but is consistent 
with the recent LHC results with 125 GeV Higgs-like boson and without SUSY particles j^. 

It would be interesting to study further constraints for R-axion with relatively large mass. 
When the R-axion mass is larger than 1 GeV, various decay channels to hadronic particle 
open. We expect that weaker but non- negligible constraints for large decay constant will be 
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imposed, thought analysis would become involved. 

A phenomenologically viable model with long-lived R-axions can be constructed by in- 
troducing a mechanism diluting the R-axion density in the early universe. Although our 
primary interest was gauge mediation models, it would be easy to apply our analysis the 
closely related situations such as spontaneous R-symmetry breaking in thermal inflation 
models [61 1. 



In this paper, we have not explicitly shown a mechanism of vacuum selection of false 
vacuum. Existence of the R-string and walls are highly depend on the scenario of the early 
stage of universe. Also, as mentioned in the Introduction, imhomogenious vacuum decay by 
impurities such as a cosmic string depends on the details of the scenario 28|. So it may be 
useful to show an explicit example of full scenario and study R-axion cosmology in detail. 
This is beyond the scope of our study, so we will leave it as a future work. 



Acknowledgment 

The authors would like to thank M. Eto for the early stage of collaboration, and T. Hi- 
ramatsu, Y. Inoue, K. Nakayama, A. Ogasahara, A. Ringwald, and T. Takahashi for useful 
comments and discussions. KK would like to thank Kyoto University for their hospitality 
where this work was at the early stage. TK is supported in part by the Grant-in-Aid for 
the Global COE Program "The Next Generation of Physics, Spun from Universality and 
Emergence" and the JSPS Grant-in-Aid for Scientiflc Research (A) No. 22244030 from the 
Ministry of Education, Culture, Sports, Science and Technology of Japan. YO's research is 
supported by The Hakubi Center for Advanced Research, Kyoto University. 



Appendix A: Higgs sector 

In this appendix, we show the Higgs sector and its mixing with the R-axion following 
Ref. Ml- Here, we concentrate on the neutral components, and H^, of the Higgs sector 
in the minimal supersymmetric standard model. The R-axion appears through the so-called 
B-term. Then, the relevant terms in their scalar potential are given by 

-(e*'^/(^^")E/xi7M + c-c.), (Al) 
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where /i is the supersymmetric mass, originated from the /i term, rnjj^ and are soft 
SUSY breaking masses squared for Hu and Hd, and Bfi is the SUSY breaking B-term. The 
B-term has an R-charge, and the axion appears there. 

At the potential minimum, the Higgs fields develop their vacuum expectation values and 
the electroweak symmetry is broken. Around the vacuum, we decompose the neutral Higgs 
fields as 



(A2) 



where and Vd are VEVs of and H^. We denote = "^u + "^d' which is related to the 
Z-boson mass mz as v'^ = 4m|/(5f^ + g'"^) = (246)^ (GeV)^. Also we denote their ratio as 



tan /3 = — . 

Furthermore, the stationary conditions, 

dV dV 



(A3) 



0, 



(A4) 



at = Vu and = Vd, lead to the following relations: 

+ = Bfi cot (3 + 2^^! cos 2/3, 

l/xp + = i?/itan/3 — -m.i cos 2/3. 

2 



(A5) 



Using them, the mixing mass matrix of the axial parts, C,u,d and the R-axion is given by 



cot (3 1 —r cos (3 

1 tan /3 — r sin /3 
.2 , 



^ — r cos (3 — r sin (3 r sin (3 cos /3 j 

Then, the mass eigenstates are given by 

/ 





Sn 
id 



(A6) 



V ^ / 



sin /3 
Kcos/3 



- cot /3 
Ksin/3 — Kr sin /3 cos /3 



Kr cos^ (3 sin /3 /tr sin^ (3 cos /3 



(A7) 



where Gq and a denote the would-be Nambu-Goldstone boson and low-energy R-axion, 
respectively. 
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Appendix B: R-axion production from R-saxion decay 

Here we estimate the R-axion abundance from the R-saxion decay and show that it is 
smaller than those from R-axion dynamics. 

Since we have assumed that the R-symmetry is restored in the early Universe, there 
should be homogeneous R-saxion oscillation associated with the spontaneous breaking of 
R-symmetry. It can take place when the Hubble parameter becomes smaller than the saxion 
mass. Here we assume that the R-symmetry is broken aX H = rris. If R-saxions receive 
thermal mass, the Hubble parameter at the time of phase transition becomes lower, but it 
requires large reheating temperature and we do not consider it here. The energy density of 
R-saxion is given by 



where the subscript "so" indicates that the parameter or variable is evaluated at the onset 
of R-saxion oscillation. The energy density of R-saxion oscillation decreases as Ps,osc oc 
due to the Hubble expansion, and gradually R-saxion decays into R-axions at H^^ = Fgax- 
Here the subscript "sd" represents the parameter or variable is evaluated at R-saxion decay. 
The number density of R-axions from saxion decay, then, is evaluated as 



where s(T) = {2ng^s(T)/4:5)T^ is the entropy density. Such R-axions are relativistic at 
R-saxion decay and loose their energy due to the cosmic expansion. After some time, 
they become nonrelativistic. The energy-to-entropy ratio, Pa.sax/s, is fixed at that time 
(if reheating is completed.) Therefore, we can estimate the abundance of R-axions from 
R-saxion decay as 



with Tso = (7r^5f*(Tso)/90)^/^(msMpi)^/^. Here we neglected the interaction of R-axion and 
assumed that the number of R-axions in a comoving volume is conserved. We can easily 




(Bl) 




(B2) 




(B3) 
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show that the abundance of R-axions from R-saxion decay is always smaller than that from 
R-string and R-string-wall system. Note that if the phase transition is driven by thermal 
potential, R-axion abundance from R-saxion decay is larger than that is estimated above. 
However, as noted, it requires high reheating temperature, in which the model has already 
been constrained by the thermal contribution strictly. As a result, the conclusion does not 
change. 
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